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1.  INTRODUCTION 

In  the  sciences  and  other  disciplines,  it  is  quite  common  to 
encounter  situations  where  random  variables  appear  in  combinations. 
Two  combinations  are  the  product  and  the  ratio  of  two  random  vari¬ 
ables  .  Given  such  combinations ,  it  is  required  to  know  the  mean  and 
the  variance  and  it  is  preferable  that  the  mean  and  variance  be  deter¬ 
mined  from  the  original  random  variables.  Thus,  if  U  =  XY  and 
V  =  X/Y  then  we  want  to  know  the  arithmetic  means  and  variances  of  U 
and  V  based  on  our  knowledge  of  X  and  Y. 

2.  ON  THE  PRODUCT  OF  TWO  RANDOM  VARIABLES 

The  mean  and  the  variance  of  a  product  are  well  known  [  5  ]  . 
To  obtain  the  result  in  a  more  simple  manner,  we  note  that 


(2.1) 


Cov(X.Y)  =  E(XY)  -  E(X)  E(Y) 


provided  that  the  moments  exist.  Henceforth,  we  will  assume  that  the 
moments  exist. 

By  rearranging  (2,1),  we  obtain 

(2.2)  E(XY)  =  Cov(X,  Y)  +  E(X)  E(Y) . 

Further,  the  definition  of  the  linear  correlation  coefficient  ^  is 
given  by 

.v  Cov(X, Y) 


(2.3) 


PX,Y 


[Var(X)l  l/2[Var(Y)]l/2 
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By  introducing  (2.3)  into  (2.2),  it  is  seen  that 

(2.4)  E(XY)  =  E(X)E(Y)  +  px  YCVar(X)]l/2[Var(Y)]l/2  . 

To  obtain  the  variance  of  a  product  of  random  variables,  we 
note  that 

(2.5)  VarOCY)  =  E(X2Y2)  -  [E(XY)1  2, 
and  that 

(2.6)  CovfX^Y2)  =E(X2Y2)  -  E(X2)E(Y2). 

Then  from  (2.1)  and  (2.5), 

Cov(X2,Y2)  =  Var(XY)  +[Cov(X,Y)  +  E(X)E(Y)1  2  -  E(X2)  E(Y2) 

=  Var  (XY)  +[Cov  (X,Y)]  2  +  2  Cov(X,Y)  E(X)E(Y) 

-  Var  (X)  Var  (Y)  -  Var(X)  [E(Y)1  2 -  Var(Y)  [E(X)1  2. 
Then ,  by  rearrangement  of  terms , 

(2.7)  Var(XY)  =  Cov(X2,Y2) -CCov(X,Y)]  2 -2  Cov  (X,Y) E(X)E(Y) 

+  Var  (X)  Var(Y)  +  Var (X)  [  E(Y)  1  2  +  Var(Y)  [  £(X)]  2 . 

If  X  and  Y  are  independent  or  uncorrelated, (2. 2)and  (2.7)  reduce  to 

(2.8)  E(XY)  =  E(X)  E(Y), 


and 

(2.9)  Var(XY)  =  Var(X)  Var(Y)  +  Var(X)  C E (Y) ]  2  +  Var(Y)  C E (X) ]  2. 

In  addition,  if  the  two  variables  are  independent  or  uncorrelated,  then 

(2.10)  Var(XY)  =  Var(X)  Var(Y) 


only  if  the  means  of  both  random  variables  are  zero. 


3.  ONE  SET  OF  RESULTS  FOR  THE  MEAN  AND  VARIANCE  OF  A  RATIO 
OF  RANDOM  VARIABLES 


(3.1) 


PY,1/X 


•3'  e| 

(i) 

-  E(Y) E 

(x) 

CVar(Y)] 

1  1/2 

P^fxl 

1/2 

Then  as  in  (2.4), 

(3.2)  e(y)  =  E(Y)  Ey)+  PY,1/*C  Var(Y)1  ^  Var("x) 1  ^ 
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provided  that  all  the  terms  exist.  In  Section  5,  it  is  shown  that 

E[-^p]  exists  provided  that  the  density  function  has  a  zero  of  order 

at  least  r  at  the  origin.  For  example.  E[ - ]  does  not  exist  for  the 

if 

Normal  Distribution. 

We  obtain  Var  (“^')*  by  replacing  X  with  •—  in  (2.71. 

(3.3)  Var(Y)  =  Cov(Y2'^2)“C  Cov(Y#*y)]  2-2Cov^Y,-^-)e(Y)e(— ) 

+  Var(Y)  Var(y)  +  Var(Y)[  E(-^-)l  2  +Var  (-^-)  [  E(Y)]  2. 

If  Y  and  ~  are  uncorrelated,  or  if  Y  and  X  are  independent 
( see  Section  5 ) ,  then 

(3.4)  E(|-)=E<Y)E(t) 
and 

(3.5)  Var(~ )=Var(Y)  Var  (-—)  +Var( Y)  e(y)2  +Var {jf)  [E(Y)1  2' 

under  the  usual  assumptions  of  existence. 

If  Y  and  X  are  uncorreiated  Normals,  E^-^-)  and  Var(*x^)  do 
not  exist.  The  proof  of  the  lack  of  existence  of  the  mean  is  shown  in 
Section  5.  Since  the  mean  does  not  exist,  neither  does  the  variance. 

In  all  of  the  above,  we  note  the  requirement  for  information 
about  the  reciprocal  of  the  denominator  term. 


4.  ANOTHER  RESULT  FOR  THE  MEAN  AND  VARIANCE  OF  A  RATIO  OF 
RANDOM  VARIABLES. 

By  using  (2.3)  we  can  write  Y 

E(Y)  -  E(X)  El—) 

S’X.Y/X'  <X' 


(4.1! 


[  Var(X)] 
Then,  by  rearrangement  of  terms. 


i/2[-  m 


1/2 


ineu,  uy  icauaiiyciiiciii  ui  uziuo,  -  /or  #  y  %-t 

.  .  [Y  \  _  E(Y)  pX,YA[Var(X>1  LVar(r)J 

LlX  /  ~  E(X)  "  E(X) 

(x,t) 


1/2 


=  im__ 

E  (X ) 


E(X) 


provided  that  the  terms  exist , 
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Y 

(Also  from  (2.7),  where  we  replace  Y  by  —  , 

(4.3)  Var(Y)  =  Cov  [x2,  (-|-)  *]-  [covfx,-^-)]  2‘ 

2  Cov  (x,  “)  E(X)  Efy)  +  Var(X)  Var  ( ■— )  + 

Var(X)  [e  (“)]  2  +  Var(“"j [  E(X)]  2. 

Rearrangement  gives  T  /  v\l9 

,  v  *  Var(Y)  -  Cov[x2,  f-)2  +  [cov(x,  J2 

(4.4)  Var  (jr)  = - - - “ - o - 

X  Var(X)  +  C  E(X)1 

2  Cov  (x,  |-)  E(X)  E  (y)  "  Var(X)  [e  (”)] 

Var(X)  +  E  E(X)  ]  2 
Introducing  (4.2)  into  (4.4),  we  get 

(4.5)  Var^)21  jvar(Y) -Cov  ^2,^j- [cov(x,  ^)J 

2  E(Y)  Cov  (x,  “)  -  VarM-  |e(Y)  -  Cov  (x, 

'  x  '  Ee(x))2L  x/ 

Y 

If  X  and  —are  uncorrelated, 

(4.6)  E  (-^-j  =  Provided  that  E(X)^0, 


rar(X)4{E(X)] 


and  a  ^ 

(4.7)  Varf^U  Var(Y)  [E(X)]  -Var(X)  [E|Y)1_  # 

X  [E(X)1  {Var(X)  +  [E(X)l  } 
provided  that  E(X)  /  0. 

Analagous  with  Section  3,  we  are  highly  dependent  on  infor¬ 
mation  about  relationships  between  the  numerator  random  variable  and 
the  ratio  of  the  random  variables. 

5.  DETERMINATION  OF  THE  MEAN  OF  A  RATIO  BY  USE  OF 
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=  Ey{tYl  EX|y(  x”lY)}  • 

If  we  consider,  as  a  special  case  the  Bivariate  Normal,  we 
know  that  the  marginal  distributions  including  the  conditional  are 
Normal  distributions.  And  for  the  Normal,  it  is  well  known  that 

E[-^-]  does  not  exist  so  that  E[  —  ]  does  not  exist. 

More  generally,  if  X  and  Y  are  independent,  then 

E[y]  =  E[Y]  E[-yl 

since  if  Y  and  X  are  independent,  then  Borel  functions  of  Y  and  X 
are  independent. 

The  phenomenon  of  the  occurence  of  the  non-existence  of 

E[— ]  is  portrayed  more  graphically  in  the  case  of  independence  in  a 

paper  by  K.  A.  David  [2,  pp.  122-123  ]  as  follows: 

Let  G(x)  be  a  cumulative  distribution  function.  Define  a 
central  moment  of  negative  order  -r  by 

(5.2)  (x)  =  x  r  dG(x)  dx. 

This  moment  exists  in  the  discrete  case  and  also  in  the  continuous 
case  provided  that  the  density  function  g(x)  has  a  zero  of  order  at 

least  r  at  the  origin  (i.e.,  lim  x  rg(x)  =  0). 

x*»0 

Consider  the  case  where  the  random  variables  are  continuous . 
Let  Y  be  independent  of  X  with  density  function  f(y);  and  if 

p.  ’  (x)  and  p. '  (y)  exist, 

-r  r 

(5.3)  P^  _(x)  p.  My)  =J*°°x  f  g(x)  dx  yr  f(y)  dy. 

-f  f  —00  —00 

or, 

(5.4)  E  [(^)r]=/r/_t1'F)r9(x)f(y)dxdy=ECYr]  E^j  . 

In  order  that  the  density  function  g(x)  have  a  zero  of  order  of 
least  r  at  the  origin,  we  are  saying,  in  effect,  that  the  random 
variable  is  either  all  positive  or  all  negative,  and  if  0  is  in  the  range 
of  the  random  variable,  then  g(x)  must  be  an  infinitesimal  at  X=0. 
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6.  RAO'  S  PROCEDURE  FOR  DETERMING  E(-£) 

Rao  [7  ]  determines  the  Expected  Value  of  the  ratio  of  Normally 
distributed  variates  in  essentially  the  following  fashion: 

2  2 

Let  X  and  Y  be  bivariate  normal  (p  ,  p  ,cr  ,  cr  ,p); 

*<  +  "'  Y-|x  X-/  7  X  y 

Now,  let  U  = - and  Z  = - —  ,  so  that  Y  =  p  +  U  a 

a  a  y  y 

y  x 

and  X  =  p  +Zor  .  From  the  definitions  of  U  and  Z,  we  have, 

X 

E(U)  =  E(Z)  =  0  and  2  =  a2  =  1. 

u  z 


Then, 


(6.1)  f  = 


p  +  U  a 

_Y _ Y 

ll  +  Z  a 

^X  X 


A  nv  /  nx  Ux  +  Z", 


p 

_  _Y_ 


/  Uff  \  /  Za  \  -1 


Rao  expanded  the  last  term  under  the  assumption  of  convergence  of  the 

sum.  It  is  well  known,  however,  that*——  converges  provided  that 

•i  +a  Z  a 

-1  <  a  <  +  1 .  This  requirement  in  our  case  is  that  -1  < - —  <  1  or 


equivalently  Zct  <  p  I.  In  terms  of  the  original  random  variable, 


the  requirement  is  that  -1  < 


x-  k 


<  +  1  or  that  0  <  X  <  2  p  for 

X 


p  >  0,  and  that  2p  <X<0  for  p  <0.  But  since  -» <X  <+  <»  for 

X  XX 


the  Normal  Distribution,  the  requirement  is  not  satisfied.  The  pro¬ 
cedure  requires  that  p^,  p  ^  0,  also. 

Nevertheless,  the  procedure  is  appropriate  for  random  variables 
which  satisfy  the  convergence  conditions  so  we  continue  the  develop¬ 
ment  of  Rao' s  procedure,  where  Z  and  U  are  appropriate  standard¬ 
ized  random  variables. 


We  note  that  finite  limits  of  the  denominator  random  variable 
are  a  prerequisite.  As  examples,  the  procedure  is  valid  for  approp¬ 
riately  parameterized  uniform  or  binomial  random  variables.  It 
probably  is  appropriate  for  cases  where  we  say  that  the  Normal 
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distribution  satisfies  the  data,  since  we  probably  mean  that  a  trun¬ 
cated  Normal  satisfies  the  data,  it  being  rare  for  an  item  or  test 
measurement  to  take  on  both  positive  and  negative  values.  By  slightly 
rewriting  (6.1)  and  performing  the  expansion,  we  obtain 

Z  a  Z*  a2  Z3  a3 


(6.2) 


X 

X 


“x  '  y  y'  l  “x 


LL 

_  _JL 


a 

+  -X 


Za 


1  - 


„2  2 
Z  a 
x 


zV 

X 


P-v 


U  *• 


ZU  or 


Z2  U  a2 
x 


Taking  Expected  Values ,  we  obtain 


(6.3) 


tl 


it 


i  + 


u3E(Z2) 


IX 

1 


ffE(ZJ) 

x 


(-1  )m  a™  E(Zm) 
m 

a2  E(Z2U) 
x  _ 


+  •  •  • 


y  - 


cr 

_X_ 


a  E(ZU) 
x _ 


= 


00 

1+  2 
1=2 


(-1)' 


(-l)1*-1  an  E(ZnU) 

X 

n 

^x 

•,  4E(zi) 


«,  (-l)j  oj  E(ZjU) 

+  ff  S  - TTT - 

y  j=i  J+i 


i 

^x 


If,  in  addition  to  the  prior  conditions  on  the  range  of  X  and 
that  ^  0,  we  specify  that  Z,  U  are  bivariate  central  symmetrically 
distributed.  Then  from  [8,pp.23],  we  have  that 


-f 
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21  +1 

E(Z  )  =  0,  i  ~  0,1, .  .. 


E(Z  U)  =  0  for  i  an  even  integer. 


Consequently, 


(6.4)  E 


fxl3 

lx  J  a 


00 

l  +  z 

j=l 


c J2 *  E(Z2j) 
x 


y<  z 

y  j=i 


a2j“1E(Z2j_1U), 

x 


If,  in  addition,  we  add  the  condition  of  independence 
between  Z  and  U  or  X  and  Y ,  then  (6.4)  reduces  to 


(6.5) 


,  pc  i .  5 

'LX  J  LI. 


00  2i 
1  +  Z  or 
i=l  x 


which  exists  under  the  usual  conditions . 

We  note  that  (6.3),  (6.4)  and  (6.5)  presented  above  could 

Y 

have  been  obtained  by  use  of  a  Taylor  Expansion  of  —  ,  the  classical 
Propagation  of  Errors  procedure . 

We  conclude  this  Section  by  noting  that  Rao’ s  procedure  for 
determining  E^~jis  valid  only  under  very  special  circumstance. 

These  are  not  satisfied  by  the  Bivariate  Normal. 

7.  THE  VARIANCE  OF  ~  ,  FOLLOWING  FROM  RAO’  S  PROCEDURE 


Following  the  classical  definition,  we  can  write 
/  \  f  2 


(7.1) 


Var(i)=  e[t-  E(i)]  • 


For  the  terms  in  the  square  brackets,  we  utilize  (6.2)  and 
(6.3)  and  after  extensive  algebraic  simplification,  we  arrive  at 

<7-21  MiY- 2-fr  +  j^-ir  Var<z‘>  + 

|Jx  l  ^x 
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(7.4) 


00  00 
2  2  2 


a 

'  4  X 

00 
i  r« 

2  )  2 
<  ( 

+  £ 
i=2 

al+k 

X 

( -1 )l 

i+k 

^x 

Var(Z  )  + 


,i  „k, 


8 


1=1  M? 


Var  ( ZJ  U ) 


/ 

)' 


and  the  approximation  reduces  to 
(7.5)  Var  (~x~)^ 


2  2,  22 
a  (J  +  a  o 
-  X  *x  y 


8.  THE  DISTRIBUTION  OF  THE  RATIO  OF  NORMAL  VARIATES 

C.  C.  Craigtl,  pp.  24-32]  presents  Fieller’ s  development 
[3,  pp.  424-440  ]  on  the  distribution  of  a  ratio  and  Geary’ s  approxi¬ 
mation  [4,  pp.  442-460  ]  .  The  following  is  a  summary  of  the  results: 

2  2 

Assume  (X,Y)  is  Bivariate  Normal  or  BVN  (p.  ,p.  ,a',a  ,  p). 
Thor.  x  y  x  y 


Then 

(8.1)  f(x,y)  = 


“  2  p 


_ I _  -1/2(1-  p*\ 

O  /I  2.1/2 

2  7r(l-  p  )  axay  I 


y  i 

Let  v  =  and  let  x  =  x.  Then,  dv  dx  =  ~r^r  dx  dy.  We  define 
x  |x( 

(8.2)  g(v)  =/+w|x|f(x,xv)  dx  =  L.  x  f(x,xv)  dx  -  x  f(x,xv)  dx 

—00  (J  —00 

=  /+°°x  f(x,xv)  dx  -  2  x  f(x,xv)  dx 
—00  —00 

Letting 

|  QQ  Q 

Q(v)  =/  x  f(x,xv)  dx,  and  R(v)  =  /  -2x  f(x,xv)  dx,  then 

'—00  -00 
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(8.3) 


g(v)  =  Q(v)  +  R(v) 


where  R(v)  >  0  since  -2  JU  x  f(xv,x)dx>0. 


First,  it  can  be  shown  that 

M* ~  p  VV(HV  +  VPX>  +  VPV<' 
(8.4)  q(v)=  y  y — | TatH 

v  2  7T  ( ff  -  2  p  v  a  a  +  v  a  )' 
y  r  x  y  x 

1  f  (|xv-V|ix)2  ) 


2  2.  ,22 
0  -2pva  a  +v  a 
L  y  H  x  y  x 


Now,  R(v)  is  of  the  same  form  as  Q(v)  with  the  exception  that  the 
limits  of  integration  are  different.  Thus 

(8.5)  R(v)  =  -2Q(v)*(^)  . 


where  $  f— r~ “)  =  /  A<T  ~~z 

\  to  i  -«  JTn 


\  to  J  ^JTT 

Now,  returning  to  the  distribution  of  the  ratio  (8.3), 
g(v)  =  Q(v)  +R(v),  we  replace  R(v)  by  the  result  in  (8.5)  to  obtain 


(8.6)  g(v)  =  q(v)  -  2  Q(v)  $ 


=  Q(v)  1-20 


f— 1 

x  On2  LA  J 

J-jl)+  «•  2J 

H  Aaj  ir  /A" 

-  —  [—1 
,  „V)  2a2  UJ 

f -ill  +  -HS - — 

l  Acr/)  7T  JT 


"'I  \tofj  nJT 

Thus  we  have  the  distribution  of  the  ratio  in  closed  form,  where  Q  (v) 
is  given  by  (8.4),  and 
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(8.7) 


2  2  2 
A  =  a  -2 pva  a  +  v  a 
y  r  x  y  x 

B=tlx‘'y  +  vVx'p(py+Vtlx)  Vy 

2  2  2  2 
C  =  a  cr  -  2  p  a  a  y.  ll  +  u.  a 
*  x  y  *  y  *  x” y  *  y  x 

2  (1-  p2)  a2  a2 

cr  =  ^  x  y  , 


we  can  rewrite  Q(v)  by  using  (8.7)  as  follows: 

AC  -  B2 1 


(8.8)  Q(v)  = 


B 


/27a  3/2 


( 


Now,  we  return  to  (8.5),  and  (8.7).  We  note  that  as  u  ->», 

X 


th"  n  B-*  oo  and 

C->*oo  .  As  a  result 

1  C 

(8.9) 

e 

20  A  .*0  , 

and 

"’(-it)-*’  • 

From  this,  we  conclude  that 

(8.10) 

R(v)-*0  as  ux-> 00  • 

Therefore,  for  u  sufficiently  large,  we  can  get  a  good 
approximation  of  g(v)  by  Q(v),  only. 

Thus,  for  the  case  where  y.  is  large,  from  (8.4)  we  have 

2  2  ~  2 
n  a  -  per  a  (u.  +  vg,  )  +vll  a 

*  ,y/?xe 

/27T((7-2pV<7  a  +  CT  )  * 
y  ^  x  y  x' 

We  impose  the  transformation 


o  _  _  i-  c 

a  -2  ova  a  -  v  a 
(.  y  r  x  y  x) 


(8.12) 


z  = 


vll  -  ix 


,  2  .  2  2.1/2  ' 
(a  -  2  pva  a  +v  a  ) 
y  r  x  y  x 


with  dz  = 


B 


,3/2 


dv, 
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and  we  note  that  the  denominator  of  Q(v)  is  positive  for  all  v. 
Hence 

vu.  -  a 

(8.13)  Q(z)=Q  I - - - X- 


,  2  .  2  2,1/2 

.  (0  -2  pva  a  +  v  a  ) 

L  y  x  y  X 

It  follows  that,  as  for  a  finite,  z  approaches  a  N(0,1). 

This  is  the  Geary  result  [4  ]  ?  X 

If  in  (8.6),  u  =  u  =0,  we  obtain 
x  y 


(8.14)  gx(v)  * 


a  -  p  2)1/z 


\  x 


2  Pv  +  v 


2  ffx 


a  generalized  Cauchy  distribution. 


The  cumulative  distribution  function  for  (8.14)  is 


(8.15)  G(v)  =  g(v)  dv  =  arc  tg 


a  v-  pa 
x  r  y 


i  ,  i  2,1/2 

IV1- p  > 


as  obtained  from  Gradstein  and  Ryzhik,  [6,  p.  82  ]  . 


The  generalized  Cauchy  obtained  by  this  writer  is  a  simple 
extension  and  has  utility  more  from  a  theoretical  viewpoint  than  an 
applications  viewpoint. 


However,  the  Geary  result  has  great  applicability  as  do  the 
approximation  results  in  Sections  6  and  7  when  one  is  willing  to 
assume  that  X  and  Y  are  Bivariate  Normal.  As  noted  in  Section  6, 
in  the  real  world,  few  if  any,  random  variables  are  truly  normally 
distributed  since  few,  if  any,  can  range  from  -  » to  +  «  .  Few,  if  any, 
will  take  on  both  positive  and  negative  values!  Thus  the  mathematical 
abstraction,  the  Normal  Distribution,  can  be  at  best  a  good  approxi¬ 
mation,  and  a  good  approximation  only  around  the  mean  of  the  data, 
since  ordinarily  one  does  not  see  data  points  in  the  extreme  tails  of 
the  Distribution. 


3^3 


FRISHMAN 


9.  SUMMARY 

In  Sections  3  and  4,  exact  expressions  are  given  for  the  Mean 
and  the  Variance  of  a  Ratio  of  random  variables ,  regardless  of  the  form 
of  the  random  variables .  One  problem  is  that  certain  correlation  co¬ 
efficients  must  be  known  or  assumed  to  be  zero  in  order  to  use  the 
formulae.  Also,  there  is  the  usual  requirement  that  the  moments  exist. 
Based  on  the  results  of  Section  5,  we  know  that  these  exact  results 
are  not  appropriate  for  the  Normal  Distribution  and  those  other  distri¬ 
butions  where  the  random  variable  takes  on  negative  and  positive 
values  and/or  the  value  0  with  probability  greater  than  0 . 

In  Section  6,  a  procedure  is  given  for  getting  the  Mean  of  a 
Ratio.  The  requirement  is  that  the  denominator  random  variable  be 
such  that  either  2^x<  X<0  when  ^  is  negative  or  that  0<X<2 p_x 
when  ux  is  positive.  Also,  we  require  that  the  Mean  of  the  numerator 
Uy/0.  In  Section  7,  we  obtain  the  Variance.  Since  the  results  in 
both  Sections  require  summing  infinite  series ,  approximations  can  be 
obtained  by  terminating  the  summation  wherever  desired. 

Finally  in  Section  8,  we  sketch  out  the  development  for  an 
approximation  of  the  distribution  of  the  ratio  if  both  random  variables 
are  Normally  Distributed,  in  the  sense  that  the  Denominator  variable 
does  not  take  on  both  positive  and  negative  values.  It  is  claimed  that 
it  is  rare  for  actual  random  variables  to  be  Normally  distributed  in  the 
mathematical  sense .  Rather,  the  Normal  Distribution  is  a  reasonable 
approximation  at  best,  and  under  these  circumstances,  the  results 
presented  in  Sections  6,7  and  8  are  appropriate  to  be  used. 
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